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Abstract. We consider a particle system in continuous time, discrete population, with spatial 
motion and nonlocal branching. The offspring's weights and their number may depend on the 
mother's weight. Our setting captures, for instance, the processes indexed by a Galton-Watson 
tree. Using a size-biased auxiliary process for the empirical measure, we determine this as- 
ymptotic behaviour. We also obtain a large population approximation as weak solution of a 
growth-fragmentation equation. Several examples illustrate our results. 
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1. Introduction 



In this work, we study the evolution of a Markov process indexed by a tree in continuous time. 
The tree can represent a population of cells, polymers or particles. On this population, we 
consider the evolution of an individual characteristic. This characteristic can represent the size, 
the age or the rate of a nutriment. During the life of an individual, its characteristic evolves 
according to an underlying Markov process. At non-homogeneous time, the individuals die 
and divide. When one divides, the chai'acteristics of the offspring depend on those and their 
number. This model was studied in ||2j|5]|6l|71[T5]|20l. Here, we study the asymptotic behaviour 
of the empirical measure which describes the population. Following Q, we begin to prove a 
many-to-one formula (or spinal decomposition, tagged fragment ...) and then deduce its long 
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time behaviour. This formula looks like the Wald formula and reduces the problem to the study 
of a "typical" individual. Closely related, we can find a limit theorem in discrete time in [11], 
in continuous time with a continuous population in |[T6ll and for a space-structured population 
model in |[T5l . Our approach is closer to ^ and extends their law of large number to a variable 
rate of division. This extension is essential in application [6|. In our model, the population is 
discrete. It is the microscopic version of some deterministic equations studied in |[27l[36l[37ll . 
Following |[T7ll43l . we scale our empirical measure and prove that these P.D.E. are macroscopic 
versions of our model. Before expressing our main results, we begin by giving some notations. 
If we start with one individual then we will use the Ulam-Harris-Neveu notation |i5j| : 

• the first individual is labelled by 0; 

• when the individual u divides, then his K descendants are labelled by ul, uK; 

• we denote by T the random set of individuals which are dead, alive or will be alive; 

• it is a subset ofU = Um>o (N*)™, where N = {0, 1, ...} and (N*)0 = {0}; 

• we denote by Vt the set of individuals which are alive at time t; 

• for each u ^ T, a{u) and (3{u) denote respectively the birth and the death date of the 
individual n; 

• we denote by Nt the number of individuals alive at time t; 

• for each u and t S [a{u), (5{u)), the characteristic of the individual u is denoted by 

Xf. 

The dynamics of our model is then as follows. 

• The characteristic of the first individual, (^t')tg[o,/3(0)) is distributed according to an un- 
derlying cadlag strong Markov process {Xt)t>o- For sake of simplicity, we will assume 
that X = {Xt)t>Q takes values in a closed subset E of and is generated by 

Gf{x)=h{x)-Vf{x) + aAf{x), (1) 

where o? G N*, 6 : M*^ — >■ M'^ is a smooth function and a G M+. Here, / belongs to the 
domain ^{G) of G. This domain is described in Section |2] Note that our approach is 
available for another underlying Maitov process. 

• The death time /3(0) of the first individual verifies 

¥[m>t\xt,s<t) = j\{xi)ds, 

where r is a non negative, measurable and locally bounded function. Notice that a(0) = 
0. 

• At time /3(0), the first individual splits into a random number of children given by an 
independent random variable K of law {pk)k£N*- We have a(0) = ... = a{K — 1) = 

my 

• We assume that the mean offspring number, which is defined by m : x J2k>o ^Pk{x), 
is locally bounded on E. 

• The characteristics of the new individuals are given by {Fj^^ (-'^^(0)_ ) ®))i<j<K, where 

is a uniform variable on [0, 1]. The sequence (-Fj'^^)j<A:,fc6N* is supposed to be a 
family of measurable functions. 

• Finally, the children evolve independently from each other like the first individual. 
The last point is the branching property. To obtain a limit theorem, we follow the approach of 
|[5l . In this paper, the cell's death rate r and the law of the number of descendants {pk)k>i are 
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constant. A many-to-one formula is proved: 

1 



dO. 



(2) 



(3) 



ueVt 

where Y is generated, for any smooth function / and x E -E, by 

Aof{x) = Gf{x)+rmY,^ [' ( i ^ 0)) - /(x) | 

This process evolves as X, until it jumps, at an exponential time with mean 1/rm. We observe 
that r is not the jump rate of the auxiliary process. There is a biased phenomenon. It is described 
in milOl and their references. We can interpret it by the fact that the faster the cells divide, the 
more descendants they have. That is why a uniformly chosen individual has an accelerated rate 
of division. A possible generalisation of dJ]) is a Feynman-Kac formula as in ll20l : 



E 



E 



ueVt 

where Y is an auxiliary process starting from xq and generated by ([3]). Using a Poisson point 
process, in ||6l, we get also another representation of the empirical measure to prove the ex- 
tinction of a parasite population. However, it is difficult to exploit these formulas. Inspired by 
|[T6llT7l[36l[37l . we follow an alternative approach. In Q, Y can be understood as a uniformly 
chosen individual. The problem is: if r is not constant then a uniformly chosen individual does 
not follow Markovian dynamics. Our solution is to choose this individual with an appropriate 
weight. This weight is the eigenvector V of the following operator: 

EE /' f{F^'\x,9))d9pk{x)] - fix) 

It is not the generator of a Markov process on E. It is described in the next section. Under some 
assumptions, we are able to prove the following weighted many-to-one formula: 



gf{x)=Gf{x)+r{x) 



nEuev.ViXf 



,E J2fi^t)vixr) =nfiyt)], (4) 

_uGVt 

where Y is an auxiliary Markov process, starting from xq- It is generated by ^ = M + J, where 
M describes the motion between the jumps and is defined by 



Mfix) 



G{fxV){x)-fix)GV{x) 
V{x) 



Gfix)+a 



VV{x).Vf{x) 
V{x) 



and J describes the jump dynamics and is given by 

J fix) = A{x 





f 


[Fj'\x,e)) 


dOpkix) 


EkeNEUloy{Fj'\x,e)) 


dOpkix) 



fix) 



where 



A(rr) 



Y.j:fy{Frix.e)) dO 



Pkix) 



r{x) 

n^) 
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These formulas seem to be complicated but they are very simple when applied. We also observe 
a biased phenomenon. But contrary to the previous formulas, the bias is present in the motion 
and the branching mechanism. This bias has been already observed in another context [16J. 
Also note that we do not assume that Aq is the first eigenvalue. So, it is possible to have different 
many-to-one formulas as can be seen in Remark 1377] We can find some criteria for existence of 
eigenelements in lH [T2j |34j [38l and theirs references. If Y is ergodic with invariant measure vr 
then Formula dU gives 



lim 



-E Y.fiXt)V{Xf) 

u&Vt 

for all continuous and bounded function /. We improve this result: 



Theorem 1.1 (Long time behaviour of the empirical measure). If the following assumptions 
holds, 

• Xq is deterministic; 

• the system is non explosive; namely Nt < +oo a.s. for all t > 0; 

• there exists {V, Aq) such that QV = XqV and V > 0; 

• Y is ergodic with invariant measure it; 

then for any measurable function g such that: 

• there exists C > 0, such that for all x £ E, \g{x)\ < CV{x); 

• there exists a < Aq, such that K[V'^{Yt)] < Ce"* and 



E 



V{Yt) 
then we have 



I' E E Pk{Ys)v{F!^'\Y„e))v{Fi'\Y„0))de 



< Ce 



at. 



lim e-^«* y g{X^) = W f ^dTr, 

where W = limt^+oo e-^"V(2;o)"^ J2ueVt ^i^t) and the convergences hold in probability. 
If furthermore V is lower bounded by a positive constant then 

t^^oo E 9{Xt) = Ivi/^o / y ^^// in probability . 



If r and {pkjk&i are constant then F = 1 is an eigenvector, and so this theorem generalises 
m Theorem 1.1]. On the other hand, our model is microscopic and is a scaled version of some 
deterministic models. More precisely, let {7it)t>Q be the empirical measure. It is defined, for all 
i > 0, by 

u&Vt 

Now, let Z'^") be distributed as Z and let us consider the following scaling X^") = ^'ZS'^\ We 
have: 



Theorem 1.2 (Law of large number for the large population). If the following assumptions hold 
• T > 0; 
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r is upper bounded; 

there exist A; > such that = 0/or all k > k; 

either E C M. and F^^\x, 9) < xfor all j < k and 9 G [0, 1] or E is compact; 

The starting distribution Xq*^^ converges in distribution to Xq G M. {E), embedded with 
the weak topology; 
we have 



supE [x(")(^) 



< +00. 



n>0 

then X^"^ converges in distribution in B([0, T], AA{E)) to X which verifies 

f f{x)Xt{dx)= f f{x)Xo{dx)+ f f gf{x)Xs{dx)ds. (5) 
Je Je Jo Je 

Here, B([0, T],A4{E)) is the space of cad-lag functions embedded with the skohorod topology 
ll9j|22|. We observe that if Xq is deterministic then X^ is deterministic for any time t > 0. The 
equation ^ can be written as 

k 

dtn{t, x) + V {b{x)n{t, x)) + r{x)n{t, x) = adxxn{t, x) + '^'^ Kj {r x pi. x nit, •)) . 

k>lj=l 

where X^ = n{t,x)dx and K!- is the adjoint operator of / i— /q^ f{F^^\x,9)d9. This equa- 
tion was studied in ETl [36l 1371 and Theorem II .H is relatively close to their limit theorems. We 
will see in the next section that it is also the Kolmogorov equation associated to Z. So, we 
observe that X is equal to the mean measure of Z; that is / H> E[/£; f{x) Zt{dx)]. This average 
phenomenon comes from the branching property. After a branching event, each cell evolves in- 
dependently from each other, there is not interaction or mutation. Another reason is the linearity 
of the operator Q. From the many-to-one formula, we also deduce that, in large population, the 
empirical measure behaves as the auxiliary process. The proof is based on the Aldous-Rebolledo 
criterion ||24]|42] and is inspired by OH] 13. 

In the end of the paper, these two theorems are applied to some structured population models. 
Our main example is a size-structured population. In this example, the size of cells grows lin- 
eaiiy and if a cell dies then it divides to two descendants. Thus, there is motion between the 
branching events and discontinuity during division. This model is a branching version of the 
well known TCP windows size process llT0l[T9ll28l[35l . For this example, we are able to give 
some explicit formulas of the invariant distribution, the moments or the rate of convergence. We 
also prove that, in large population, the empirical measure behaves according to the detenninis- 
tic equation ([Sjl plus a Gaussian noise. 

Outline. In the next section, we introduce some properties of the empirical measure. In Section 
[3j we focus our interest on the long time behaviour. We prove some many-to-one formulas and 
deduce a general limit theorem which implies Theorem 11.11 Section |4] is devoted to the study 
of large populations. In this one, we prove Theorem 11.21 Note that Section |3] and Section |4] are 
independent. In Section [51 we give our main example, which describes the cell mitosis. More- 
over, we give two theorems for the long time behaviour of our empirical measure in addition to 
some explicit formulas. We also give a central limit theorem for asymmetric cell division for 
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the macroscopic limit. In section |6j we finish by two classical examples which are branching 
diffusions and self-similar fragmentation. 

2. Preliminaries 

In this section, we describe a little more the empirical measure {Zt)t>o- We recall that 

yt>0,Zt=Y ^x^. 

u&Vt 

It belongs to in the space D(M+, A^(£')) of cad-lag functions with values in M{E), which is 
the set of finite measures on E. Let us add the following notations: 



ZtU) = f f{x)Zt{dx) = Y 



and 

Zt(l + xP) = / 1 + Zt(dx) = E 1 + ^^tY- 

We can describe the dynamics of the population with a stochastic differential equation. That is 
given, for any smooth function / : (x, i) i— )■ /(x, t) = /j(x) on £' x M+, by 

ft 



Zt(/j) =Zo(/o) + / / Qfs{x) + dt!six)Z,{dx)ds 
Jo Je 

^ I L ^ r i^{y.eVs-,i<r(x^_)} 



JWxK+xN*x[0,l] 
k 

(fe) 
j 



{Yfs{Fj'\x^^,9)) - fs{X^^))] p{ds,du,dl,dk,d9), 



where {B'^)u^u is a family of independent standard Brownian motions and p{ds, du, dl, dk, dO) 
is Poisson point process on M+ xU x M+ x N* x [0, 1] of intensity 

p{ds, du, dl, dk, dO) = ds n{du) dl dpk dO. 

It is also independent from the Brownian motions. We have denoted by n{du) the counting 
measure on U and ds, dl, dO are Lebesgue measures. A necessary and sufficient condition for 
the existence of our process is the non-explosion of Z: 

Assumption 2.1 (Non explosion). For all t > 0, Nt < +oo a.s.. 

For instance, we have 

Lemma 2.2 (Sufficient condition to non explosion). If r < r and pk = 0/or all k > k, where 
f,k>0, then Assumption \2n\ holds. Moreover, for any T > 0, we have 

\ft < T, E[Nt] < E[No] ^(^-1)^^. 

Proof. In this case, we can bound Nt by a branching process independent of the underlying 
dynamics. □ 
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If there is no explosion then we have a semi-martingale decomposition. Before giving it, we 
give some precisions and notations about the domain V{G) of G. If G is given by ([T]) then we 
have T>{G) = Cl{E, M), which is the set of bounded and functions with bounded derivatives. 
We denote by 'Dt{G) the domain of the Markov process {Xt,t)t>o . Also, we have Vt{G) = 
C^' (M4. X E, M). It is the set of bounded functions, which are in their first coordinate, in 
their second, and which have bounded derivatives. Finally we denote by 'D(G^) (resp. Vt{G'^) 
) the set of function / G V{G) (resp. / G VtiG) ) such that f G P(G)(resp. / G VtiG^) ) . 

Lemma 2.3 (Semi-martingale Decomposition). If Assumption 12. 1 1 holds, then for all bounded 

f = ift)t>o £ T^tiG) and t > 0, we have 

Ztift) = Zoifo)+Mtif) + Vt{f) 

where 

Ytif) = f'zs{Gfs + dsfs)ds, 
Jo 

and if f G T){G'^) then the bracket of'M.t{f) is given by 

(M(/))i = f'G{f^){x) - 2Mx)Gfs{x)Z,{dx) 
Jo 

+ I r{x) r |E/^(^r(^'^))-/^(^) I Pk{x)d9Zs{dx)ds 
Je Jo ^gpj, \^-^^ J 

Proof. It is an application of Dynkin and Ito formulas, see for instance |[23l Lemma 3.68 p .487] 
and 1221 Theorem 5. 1 , p.67] . □ 

We define the mean measure (zj)t>o, for any continuous and bounded function / on E, by 

Vt>0, Zi(/)=E(Zt(/))=E ^/(Xr) 

ueVt 

Corollary 2.4 (Evolution equation for the mean measure). Under Assumption \2.1\ if f is a 
continuous and bounded function on E then for any t > 0, we have: 

zt(/)=zo(/)+ f^s{Gf)+ I r{x) r f(Fl''\x,0)) d9 pk{x) - f{x) z,{dx)ds. 

Jo Je fc>ij=i-^o ^ ^ 

Furthermore, (zt)t>o is the unique solution of this integro-differential equation. 
The previous equation can be written as 

k 

dtn{t, x) + V {b{x)n{t, x)) + r{x)n{t, x) = adxxn{t, + ^ ^ Kj {r x pk x n{t, •)) ■ 

fc>ii=i 

where zt = n{t, x)dx and is the adjoint of / H> f^ f{Fj''\x, 9)d9. 
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Proof. We just have to prove the uniqueness. Consider two probabiUty measure-valued pro- 
cesses {fJ.t)t>o and {i't)t>o solution of this P.D.E. with same starting distribution /xq = i^o- Let 
if be a compact subset of E and let us consider the following norm 

||"2l-m2||= sup |mi(/) - ??T,2(/)|, 

where / G C^(E',M) belongs to J^k if and only if ||G/||oo < 1 and /(x) = when x ^ K. 
Now, let / G Fk, it X ^ K then Gf{x) = 0. The functions r and m are bounded on K and thus 
there exists Ck > such that ||^/||oo < Ck and then 

lAii(/) - Mf)\ <Ck[ WfJ-s - i^s\\ds. 
Jo 

Taking the supremum and using Gronwall Lemma we deduce that = ft(/) for any 

/ G J^K- As there is no explosion, we can approach any continuous and bounded function 
by dominated convergence, and then the equaUty holds for any continuous and bounded func- 
tion. □ 



3. Long TIME'S BEHAVIOUR 

Let us recall that 



gfix)=Gfix) + r{x) 



EE r f{F^'\x,e))d9pk{x)-f{x) 

k>0 1=1 



fc>0 j= 

In the following, we will prove some formulas which characterise the mean behaviour of our 
model. Then we will use them to prove our limit theorems. 

3.1. Eigenelements and auxiliary process. As said in introduction, the existence of eigenele- 
ments is fundamental in our approach. Henceforth, we assume the following. 

Assumption 3.1 (Existence of eigenelements). Assumption 12. 1 1 holds, and there exist Aq > 
and V G V{G) such that y > and QV = \qV. 

Under this assumption, we introduce the martingale {7jt{V)e~^°^)t>Q which plays an important 
role in the proof of theorem 11.11 

Lemma 3.2 (Martingale properties). If Assumption \3. 1 1 holds and 

zo(^) < +00, 

then the process (Zt(y)e~^°*)f>o is a martingale. Moreover, it converges almost surely to a 
random variable W. 

Proof. First, by corollary 12.41 we have 



zt(y) = zo(y)+ [\s{GV)ds 
Jo 

= zo{V) + Xo [ Zs{V)ds. 
Jo 



Hence for all t > 0, we have zt{V) = zqc^o*. Then if J^t = cr{Zs \ s < t} then the Markov 
properties, applied on Z, gives 

E[Zt+s{V)\J^s] = E[Zt{V)\Zo = Zs], 
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where Z is distributed as Z. Then E[Zt+s{V)\J='s] = Zs{V)e^ot ^nd thus 

E[Zt+,(F)e-^o(*+^) I J^s] = Z,(y)e^«^ 

Since {Zt{V)e~''^°^)t>o is a positive maitingale, it converges almost surely. 

Lemma 3.3 (Weighted many-to-one formula). Under Assumption \3.1\ i/Zq 
E, then we have 

= nf{yt,t)\Y^ = x^], (6) 

for any non negative function f on E x M4. and t >0, where Y is a Markov process generated 
by A. 

Proof If 7t : / zt{f x y)e^^«*y(xo)^^ then for alH > and smooth / we have 

datif) = ztiGivf) + vdtf - fgv)e-^'>'v{xo)-^ = jM/ + dtf). 

Now, by Dynkin formula, the right hand side of ([6]) verifies the same equation. The uniqueness 
can be proved as in the coroUaiy 12.41 □ 

Remark 3.4 (The first characteristic can be random). If Zq = 6^9, where is random and 
distributed according to a probability measure fi, then ^ holds, where Y starts form Yq which 
is distributed according to fx. 

Remark 3.5 (Schrodinger operator and /i— transform). Q is not a Markov generator. As we have, 
for all f smooth enough, 

gf = Bf-r{m-l)f, 
where B is a Markov generator and Q a so-called a Schrodinger operator. Its study is con- 
nected to the Feynman-Kac formula. The key point of our weighted many-to-one formula is an 
h— transform (Girsanov type transformation) of the Feynman-Kac semigroup as in |!39l|. This 
transformation is usual in the superprocesses study II16I . 

Remark 3.6 (Galton- Watson tree and Malthus parameter). If r and p are constant, then V = 
1 is an eigenvector with respect to the eigenvalue Aq = r{m — 1), where m = J2k>o^Pk 
denotes the mean offspring number. So, Zt{V) = Nt and the population grows exponentially. 
This result is already know for Nt. It is a continuous branching process |l3l|5l. Furthermore, 
since Thomas Malthus (1766-1834) has introduced the following simple model to describe the 
population evolution: 

dtNt = birth - death = bNt - dNt = XoNt ^ Nt = e^°\ 

in biology and genetic population study, Aq is sometimes called the Malthus parameter. 

Remark 3.7 (Many eigenelements are possible!). In the previous lemmas, Aq was not required to 
be the first eigenvalue. So, it is possible to have different eigenelements and auxiliary processes. 
Consider the example of |[6l, where some eigenelements are explicit; that is: 

Vx > 0, Gf{x) = ax fix) + bxf"{x), 

where a, b are two non-negative numbers. We also consider that p2 = 2 and for all j £ {1, 2}, 

E[f{F^^\x,e))]=E[fiHx)], 



□ 

= 5x0, where xq G 



rE 
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where H is a symmetric random variable on [0, 1] i.e. H = 1 — H. This example models cell 
division with parasite infection. In this case, 

gf{x) = ax fix) + bxf'ix) + r{x) {2E[f{Hx)] - f{x)) , 

where a is an eigenvalue of Q and V[x) = x is its eigenvector So, we should have 



E 



E[f{Yt)]e'''xo 



where Y is a Markov process, generated by 

Gyfix) = {ax + 26) f'{x) + bxf"(x) + r{x) {2E[Hf{Hx)] - f{x)) . 

We can see a bias in the drift terms and jumps mechanism which is not observed in ||5l|6l. When 
r is afftne, we obtain a second formula. Indeed, if r{x) = cx + d, with c > and d > a (or 
d > and c = 0) then Vi{x) = x{c/ (d — a)) + 1 is an eigenvector with respect to the eigenvalue 
Xi = d Xi > Xq = a). Thus, we should also write 



E 



-dt 



E 



fjUt) 
rUt + l 



{tXq + 1), 



where r = and U is generated by 
2bxT 



Guf{x) 



ax + 



TX + 1 



, , r(x)(Tx + 2) ( 2E[{THx + l)f{Hx)] 
f {x)+bxf (x)H — — f{x) 



TX + 1 



TX 



3.2. Many-to-one formulas. In order to compute our limit theorem, we need to control the 
second moment. As in 151, we begin by describing the population over the whole tree. Then we 
give a many-to-one formula for forks. Let T be the random set representing cells that have lived 
at a certain moment. It is defined by 

r={ueU\3t>o,ueVt}. 

Lemmas 1X8] and [X9l that follow, ai'e respectively the generalisation of IS proposition 3.5] and 
15] proposition 3.9]. 

Lemma 3.8 (Many-to-one formula over the whole tree). Under Assumption \3.I\ if Zq = 6xq, 
where xq G E, then for any non-negative measurable function f : E x ]R_|. — )■ M, we have 



E 



E/tan)-,/?W 



.ueT 



E 



fiYs,s) 



r{Ys) 
V{Ys, 



V{xo)e^'"ds 



Proof. First we have, for all u ^U, 



E 



E 



l{ner} 1^^^^ f{X^,s)r{X^)ds 

a{u) 
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because 



Thus, 



E 

=E 

=E 

=E 
=E 



{ueT} / f{X^.s)r{X^)ds 

J a{u) 

f{X^, s)r{X^)ds r{X^)e~ -f^^-^ dr 
(X;')e-/.V)'^(^")'^* dr f{X^,s)r{X^)ds 

I^(^)'^''^^'''f{X^,s)riX^)d. 
l{«er}/(^^H-,/?(«) 



{u) 

+ 00 PT 



Ja{u) 
+00 r+00 



L{«er} 



a(u) Js 
r+00 

l{«gr} / e 

Ja(u) 



E 



{ueT}f 



E 



l{ueV.}f{X^,s)r{X^)ds 



and finally, 



E 



+00 



E 



+00 



E 



Y: f{X^,s)r{X^ 



f{Ys,s) 



ds 

V{xo)e^°'ds. 



If we set g{x, s) = f{x, s) /V{x) then we have: 



E 



.«er 



+00 



□ 



E[<7(n, s)r(y,)] xE[Zs{V)] ds. 



This equality means that adding the contributions over all the individuals corresponds to inte- 
grating the contribution of the auxiliary process over the average number of living individuals at 
time s. Let {Pt)t>o be the semigroup of the auxiliary process; it is defined, for any continuous 
and bounded /, by 

Ptf{x) = E[f{Yt) \Yo = x] 

Lemma 3.9 (Many-to-one fomiula for forks). Under Assumption li.il j/Zq = 6x0, where xq S 
E, then for all non-negative and measurable function f, g on E, we have 



E f{Xr)V{X]^)g{Xl)V{Xl) 



E 



J2iVPt-J,VPt-sg)iYs) 



r{Ys) 

y{Ys) 



ds 
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where J2 is defined by 



j,{f,g){x) = J y: E Pk{x)f{F^'Hx,e))g{Fi'\x,e))de. 

^ a^b fe>max(a,6) 

The operator J2 describes the starting positions of two sibhngs picked at random. 

Proof. Let u,v &Vt be such that u ^ v, then there exists {w, u, v) G and a, 6 G N*, a ^ b 
such that u = wau and v = wbv. The cell w is sometimes called the most recent common 
ancestor (MRCA). We have 

E J2 f{Xr)ViX;^)g{X^)V{Xl') 

where u = wau and v = whv . Let Ft = (j{7is \ s < i}. By the conditional independence 
between descendants, we have 



E 



fiXr)V{X^)giX^)V{Xn 



= ^^E E lu.a^>.^f{XT)V(X]')\FR,.,^ E 
Therefore, as fi{w) is a stopping time, then using the strong Markov property and we have 



E '^{v&t}9{Xt)V{X^)\Fi3(^u,) 



E f{xr)v{xr)gix^)vm 

E E^[-'-{«'a.«'''Gr,i>/3(w)}^t-/3(«))/(-'^^('i,))^(-'^^(^)) 



uiGW a^b 



=E 



E l{^>/3(^«)}'^2(V'P^_/3(^„)/, yP^_/3(^„)S')(XJ(^)_) 6 



2\o(t-l3(w)) 



.weT 



Jo 



'j,{VPt.sf,VPt.sg){Y,)^^ 



-Xqs 



ds. 



□ 



3.3. Proof of Theorem II. 1[ In this section, we give the main limit theorem which implies 
Theorem 11.11 

Theorem 3.10 (General Condition for the convergence of the empirical measure). Under As- 
sumption \3.1\ if f is a measurable function defined on E and fi a probability measure such that 
there exists a probability measure vr, two constants a < \q and C > 0, and a measurable 
function h such that 
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(1) 7r(|/|) < +00 and'ix G E,\init^+^Ptf{x) = 7r(/), 

(2) ^x{V) < +00 and fiPtif x V) < Ce"*, 

(3) Pt\f\ <handf,Ps {J2iVh,Vh) f ) < C7e°*, 
and Zq = d^iD, where ~ /i, then we have 
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1 



lim , , 
i^+oo E[Zt{V)] 



Y,fiX;^)V{Xr) = Wxn{f), 



ueVt 



where the convergence holds in probability. If furthermore (Z4(y)e ^°*)t>o is bounded in 
then the convergence holds in L^. 

Note that the constants and vr may depend on / and fi\ Also note that Aq is not assumed to be 
the largest eigenvalue. 

Proof. As in ||5] Theorem 4.2], we first prove the convergence for / such that 7r(/) = 0. We 
have E[Zt{V)] = /i(y)e^«* and so 



E 



1 



E[Zt{V)] 



J2f{xnv{x^) 



ueVt 



E 



Ztif X y)2e-2Aot^(y)-^ =At + B-, 



where 



At = e-2^oV(^)"^lE 



ueVt 



e-^°'liiVy'E \f\Yt)ViYt) 



and 



Jo 



f{xr)v{xf)f{xi)v{x^t) 



J2{VPt^sf, VPt-sfWs, 



From (|2]i, we get limt_>+oo At = 0. Since vr(/) = 0, from ([T]), we get limt_> Ptf = 0. Then, by 
^ and Lebesgue's theorem, we obtain that, for all s > and x ^ E, 

lim J2{VPt-sf,VPt-sf){x) = Q. 

t— 5-+00 

And again by ([3]) and Lebesgue's theorem, we obtain that limt^+oo Bt = 0. Now, if 7r(/) / 
then we have 

Zt(/y)e-^°* - W^U) = ((/ - vr(/))y) e"^"* + 7r(/) (Zt(y)e-^« - w) . 

Then, as a consequence of to the first part of the proof, the first term of the sum, in the right hand 
side, converges to in L^. Moreover, the second term converges to in probability thanks to 
lemma [3^ □ 
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Proof of Theorem [77/] If / = g/V then it is a continuous and bounded function. \fh = 1 then 
all assumptions of the previous theorem hold and we get the first convergence. Now if V is lower 
bounded, we can use the same argument with g = I and / = l/V which is also a continuous 
and bounded function. □ 



4. Macroscopic approximation 

To prove Theorem 11.21 we need to use different topologies on M{E). Let {M{E), d^,) (resp. 
{M.{E),dw)) be the set of finite measure when it is embedded with the vague (resp. weak) 
topology. These topologies ai^e defined as follow. 

lim d,{Xn,X^) = ^ V/ G Co, lim E[/(X„)] = E[/(Xoo)], 

lim d^iXn^X^) = ^ V/ G a, lim E[/(X„)] = E[/(Xoo)], 

where (X„)^o is a sequence of M{E) and X^o G J\A{E). Here, Co is the set of continuous 
functions which vanish at infinity, and is the set of continuous and bounded functions. Let 
©([O, T], £■) and C{[^,T],E) be respectively the set of cad-lag functions embedded with the 
Skohorod topology and continuous functions embedded with the uniform topology [;9 |. 

4.1. Proof of TheoremlLll Let(Z(")) „>i be a sequence of random measure- valued distributed 
as Z. In this section, we consider the following scaling: X^"^ = ^Z^"^ and we describe the 
behaviour of this scaled process when n goes to infinity. 

To understand the behaviour of our model in a large population, we can consider that it starts 
from a deterministic probability measure Xo, and approach it by the interesting sequence defined 
by 



where (lfc)fe>i is a sequence of i.i.d. random variable distributed according to Xo. In other 
words, we set 

n 
fe=0 

Does X^"^ converge? Yes, it converges. Indeed, by the branching property, we have Z^") = 
X]fc=o Z^*-", where Z^*"' are i.i.d., distributed as Z and starting from Z^* = Sy^- Henceforth, if / 
is a continuous and bounded function then the classical law of large number gives 



Vt>0, lim X^")(/)=E Zl'{f) 



a.s. 



So by corollary 12.41 it implies that X^") (pointwise) converges to the solution {iJ,t)t>o of the 
following integro-differential equation: 

Mf) = M + [%s{Gf) (7) 



•^^ fc>o j=l 
Theorem 1 1 . 2 1 gives a stronger convergence. 
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Lemma 4.1 (Semi-martingale decomposition). If Assumption \2.1\ then for all f G 'D(G^) and 
t > 0, 

x1")(/) = x(")(/) + m(")(/)+v(")(/), 



where 



^t\f)= f I Gf{x)+r{x) /'5]X^/(F,('^)(x,e))-/(x)pfcd0X(«)(dx)ds, 
-'0 Je Jo frtt^ 



fceN j=i 

and m|"^ (/) is a square-integrable and cddldg martingale. Its bracket is defined by 

(M(")(/))t = - r2XW(G/2) - 2XW(/ X Gf) 
n Jo 

+ / r{x) C J2 Pk{x)d9X<i'\dx)ds. 
JE Jo ^gpj. y^^^ J 

Proof. It is a direct consequence of Lemma 1231 Indeed, if L^") is the generator of X^") then it 
verifies 

= 5tE[F;(xW)|x(") = = aiE[F^/„(zW)|z('^) = nt,]^^^ = hFf/^M, 

where Ff{ii) = F{fi{f)), F, f are two test functions and L is the generator of Z. □ 

Remark 4.2 (Non explosion). Let us recall that, by Lemma \2.2\ if the assumptions of Theorem 
\1.2\ hold then Assumption 12. 1 \ holds: that is there is no explosion. 

Let us denote by C{U) the law of U, for any random variable U. 

Lemma 4.3. Under the assumptions of Theorem \1.2\ the sequence (£(X(")))„>i is uniformly 
tight in the space of probability measures on ©([O, T], {M{E), d^)). 

Proof. We follow the approach of lITTl . According to Il42l . it is enough to show that, for any 
continuous bounded function /, the sequence of laws of X(")(/) is tight in D([0, T],IR). To 
prove it, we will use the Aldous-ReboUedo criterion. Let be the sef of functions of class 
C°° with finite support, we set S = U { 1}, where 1 is the mapping x ^ I. We have to 
prove that, for any function / G S, we have 

(1) yt>0, (xS")(/)) _ is tight; 



'n>0 

(2) for all n G N, and e,r] > 0, there exists 5 > such that for each stopping time Sn 
bounded by T, 

limsup sup P(|vg;V„(/) - Vg;)(/)| > r?) < 

rn>+oo 0<u<S 

limsup sup P(KM(")(/))5„+.-(MW(/))5j>7?)<e. 

n-s-+oo 0<u<S 
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The first point is the tightness of the family of time-marginals (x|"^(/))„>i and the second 
point, called the Aldous condition, gives a "stochastic continuity". It looks hke the Arzela- 
Ascoli Theorem. Using Lemma IZ21 there exists C > such that 



X't\f)\>k)< 



< 



E[xi"V)] 
k 

k 



which tends to as tends to infinity. This proves the first point. Let (5 > 0, we get for all 
stopping times 5„ < r„ < {Sn + S) <T, that there exist C',Cf > such that 



E[|vg(/)-V^':^(/)|]=E 



+ 



Je Jo ^gpj ^ ^ 



<C'[||G/||oo + ||/||oo] xE[|r„-5„| 
<Cj6. 

In the other hand, there exists Cj>0 such that 
E[|(mW(/))t„-(mW(/))5J] 



1, 



-E 



n 



2X(")(G/2)-2Xi«)(/G/) 



<- 



n. 



Then, for a sufficiently small 6, the second point is verified and we conclude that (X("))„>i is 
uniformly tight in D([0,T],(7W(^), 4)). ~ □ 

Proof of Theorem \L2\ Let us denote by X a limiting process of (X(")).„>i; namely there exists 
an increasing sequece {un)n>i, on N*, such that (X(""))n>i converges to X. It is almost surely 
continuous in {J\A{E)^v) since 



sup sup |xi:!)(/)-xr(/)i<- 

t>0||/||oo<l " 



k 



(8) 



In the case where E is compact, the vague and weak topologies coincide. By Doob's inequality, 
there exists C > such that 



supE 

/ 



sup 

t<T 



< 2supE [(m(")(/))t 
/ ^ 



< 



C 



n 
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where the supremum is taken over all the function / G 'D(G^) such that ||/||oo < 1- Hence, 



But as 



lim sup E 

71— > + 00 J 



Mr(/) = xr(/)-x(")(/) 



sup 

t<T 



(9) 



we have 



t f Gf{x) + T{x) f' f - fi^) P''i^)dO^'-r\dx)ds, 

0=Xi(/)-Xo(/)- rx,(G/) 
Jo 

Since this equation has a unique solution, it ends the proof when E is compact. This approach 
fails in the non-compact case. Nevertheless, we can use the Meleard-Roelly criterion BOi We 
have to prove that X is in C([0, T], {M{E),w)) and X(")(l) converges to X(l). By X is 
continuous. To prove that X'^")(l) converges to X(l), we use the following lemmas. 

Lemma 4.4 (Approximation of indicator functions). For each fc G N, there exists ipk £ ^(C) 
such that: 



\/x e E, l[fc;+oo[(a;) < ipk{x) < l[fc_i;+oo[(a;) and3C, Gipk < Cipk-i- 
Proof. See El lemma 4.2] or Ml lemma 3.3]. 

Lemma 4.5 (Commutation of limits). Under the assumptions of Theorem U .2\ 



□ 



lim lim sup E 



t<T 



0, 



where {ipk)k>o <^^^ defined as in the previous lemma. 



The proof is postponed after. Hence, a same computation to OTI gives us the convergence in 
]D)([0, T], {M.{E), w)). Thus, each subsequence converges to the equation (|7J. The end of the 
proof follow with the same argument of the compact case. 

We can give another argument, which does not use the Meleai^d-Roelly criterion BOj . As 



sup sup \xt\f)-^t\f)\<-, 
t>o ||/|U<i n 
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X is continuous from [0,T] to {M{E),dw). Let G be aLipschitz function on (7([0, T], {M{E),dw)), 
we get, 



|E[G(X(""))] - G(X)| < E 



sup d^(x(""\Xt) 

te[o,T] ^ ' 



< E 



+ E 



sup 

te[o,T] 



(i^(x^\x^)(.x(l-^fc)) 



sup d^(xi"")(- X (l-Vfc)),Xt(- X (1-^fc))) 

tG[0,T] ^ ^ 



+ sup (Xi(- X (1 - V'fc)),Xj) 

t6[0,T] 



According to Lemma 1431 we obtain that 



lim limsup E 



fc^+oo n-s>+oo [te[0,T] 



sup 



d.(x^u^^(•x(l-v.))) 



and 



lim sup d„(Xt(- X (1 - ^fc)), Xf) = 0. 

fc->+oo fg[o,T] 



Then, we have 



(xi"")(- X (1 - ^fc)),Xt(- X (1 -^fc)) 
4(xl"")(- X (l-^fc)),Xj(. X (1-Vfc))) 



Thus, 



lim limsup E 

A:-5>+oo n-s>+oo 



sup d^ (X^"")(- X (1 - Vfc)),Xt(- X (1 - V'fc)) 

tG[0,T] ^ 



0, 



by continuity of H> — V'fc) in ][D(A^(-E), d^). And finally, 

lim G fx("")) = G(X), 

which completes the proof. 

proof of Lemma W3\ If /i"'^ = E(X^"^(V'A:)) th^n we have 



/.r'' = E[x(,")(V^fc)]+ rE[/ GMx) 
Jo Ue 

Vfc>li=l J 

<f,^'' + C [\^''-' + f,^^'^ds. 
Jo 



ds 



□ 
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Now, by Gronwall's Lemma, iteration and mono tonicity, we deduce that 

Jo 



^ ^ n,k , r^2rr n,k—l 



Jo 



nl;^ "^duds 



fc-i 

/=0 



^ \^ n,k-l^ {CiTY 



+ C2 X 



n " ki 



l>lk/2\ 

where Ci, C2 and C3 are three constants. Thus, 



II 



+ C2 X 



kl 



lim lim sup/i[^''^ = 0. 

fc-5- + 00 71-S- + 00 

And finally the following expression completes the proof, 



E 



suplXnV'fc) 

t<T 



JO 



sup|M('^)(^fc)| 



□ 



5. Main example : A size-structured population model 

Let us introduce ouif main example. It is a size-structured population model which represents 
the cell mitosis. It is described as follows: the underlying process X is deterministic and linear 
and when a cell dies, it divides in two parts. Formally and with our notations, we have 



^= [0,+cx)), Gf = f' and = 1, 



(10) 



yx e E,\f 9 e [0,1], F[^\x,e) = F-'^{e)x and F^^\x , 6) = {I - F'^ {e))x , (11) 

where F is the cumulative distribution function of the random variable in [0, 1]. It verifies 
F{x) = 1 — F(l — x). In this case, one cell lineage is generated by: 



V/ G C\ Vx > 0, L/ = f\x) + r{x) [E[/(/7x)] - f{x)] , 

where H is distributed according to F. This process is sometimes called the TCP (Transmis- 
sion Control Protocol) process in computer science lITO] [T9] l28l l35l . Firstly, we prove the non 
explosion even if r is not bounded. 

Lemma 5.1 (Non explosion). Letp > 1. If for all x G R*_^, r{x) < Co{1+xP), and zo{1+xP) < 
+00, then our process is no, explosive. Moreover 



E 



sup Zs(l 

5G[0,T] 



<zo(l + xP)e^^^, 



where Cp is constant and T > 0. 
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Proof. We have 

Zt(/) = Zo(/) + /* / f{x) Zsidx) ds 
Jo Je 



+ / / Mu^v..,i<rix^ )}fiOX:_) + /((I - 9)Xf_) - fiXl^_)) p{ds, du, dl, dO) 

Using the same argument to lITTl Theorem 3. 1], we introduce r„ = inf{ t > | (l+x^) > n }; 
and we have 

tAT„ 



sup Z„(l + <Zo(l + ) + / Zs{pxP-^)ds 







t/\Tn 



+ / / l?ie\4-,Kr(X" ) 

10 JUxR+x[0,l] "~ 

(1 + + (1 - e)p - i){x^^y) p{ds, du, di, do) 

rtAT„ 

<Zo{l + xP)+ / sup Zu{l + xP)ds 

JO U&[0,sATn] 

+ 1 I ^{u€Vs-,i<r{x^_)} Pids,du,dl,de), 



lo JUxR+x[0,l] 

because (6'P + (1 - 9Y - 1) < 0. Thus there exist C > such that 

rt 



E 



sup Z„(l + x") 

we[o,MT„] 



<zo(l + xP)+ / CE 

^0 



sup Zu[l + x^) 

mG[0,sAt,i] 



ds. 



Finally, the Gronwall Lemma implies the existence of Cp such that 



E 



sup Z,(l + xP) 

se[0,iAT„] 



<zo(l + 2;P)e^f*. 



We deduce that r„ tends almost surely to infinity and that there is non explosion. 



□ 



5.1. Equal mitosis : long time behaviour. In this subsection, we establish the long time be- 
haviour of Z. We assume that 

vx > 0, V0 G [0, 1] , ) {x,e) = ix,e) = ^. 

That is, the cells divide in two equal parts. In order to give a many-to-one formula, we recall a 
theorem of ll37l : 



Theorem 5.2 (Sufficient condition for the existence of eigenelement). If there exist r , r > 
such that 

r<r<r, 

r is continuous and r{x) is constant equal to r^o for x large enough, then there exists V G 
C^(M+) such that 

c{l + x^) <V{x) <C{l + x^), 
where C, c are two constant and 2^ — 



Ao+r-oo ■ 
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So, we get a many-to-one formula with an auxiliary process generated by 

Af{x) = fix) + r{x)^^ (/(x/2) - fix)) . 
Our main result gives the two following limit theorems. 

Corollary 5.3 (Convergence of the empirical measure for a mitosis model ). If there exist r,r > 
such that 

r <. r ^ f, 

r is continuous and r{x) is constant equal to rfor x large enough, then there exists a probability 
measure vr such that, for any continuous and bounded function g, we have 



.1™ 'W Yl ai^t) = [ g d-K in probability . 



ueVt 

In particular for a constant rate r, ir has Lebesgue density: 



-oo 




x^ . e-" ™. (12) 

This explicit formula (fT2l) is not new |[36l [371 . but here, the empirical measure convergences 
in probability, while in the mentioned papers, the mean measure or the macroscopic process 
converges (see Theorem [L2] i. 

Proof of corollary \5. 3\ By Theorem 15.21 the mapping x ^ V{x/2)/V{x) is upper and lower 
bounded. Thus, the auxiliary process is ergodic and admits a unique invariant law, as can be 
checked using a suitable Foster-Lyapunov function ll33l Theorem 6. 1] (for instance, we can use 
X ^ I + x). See also lITSl . Finally, we use Theorem II. II to conclude. The explicit formula is an 
application of |f35l. □ 



We can see that the assumptions of Theorem l5.2l are strong, and not necessary: 
Corollary 5.4 (Convergence of the empirical measure when r is affine). If 

Vx > 0, r(x) = ax + b, 
where a,b >0 and a orb is positive then there exists a measure vr such that 



lim 

The convergence holds in probability and for any continuous function g on E such that Vx G 
E, \g{x)\ <C(l + x). 

Proof If r(x) = ax + b then V{x) = x ^^'^'\^°-~ - + 1 is an eigenvector and ^^2 ^4^ ^ 
corresponding eigenvalue. Henceforth, this result is a direct application of Theorem II . 1 1 □ 



Remark 5.5 (Malthus parameter). We also deduce that 

lim Nt e-^o* = W I l-dn, 

JE V 

where \q = ^^,2^4^ h Malthus parameter (see Remark WM . 
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Remark 5.6 (Estimation of r for the Escherichia coli cell). We can find some estimates of the 
division rate in the literature. An inverse problem was developed and applied with experimental 
data in ^^(see also II26II ). More recently, II13II gives a nonparametric estimation of the division 
rate. 

5.2. Homogeneous case: moment and rate of convergence. When r is constant, the process 
is easier to study since the auxiliaiy process has already been studied |[T0ll28l[35l . Here, we give 
the moments and a first approach to estimate the rate of convergence. 

Lemma 5.7 (Moments of the empirical measure). For all m G N, and i > 0, we have 
E[Zt(a;-)]=E _ 



+ 00 



ml 



X 



where Oi = 2r (1 — 2^*). In particular, 



: + -'E Eir n 



i=l \k=0 j=k,jf^i 1 



zo{dx), 



K[Zt{x)]=E 



E^r 



+00 I 

r 



X ] e ^ ZQ{dx), 



and 



E[Zt(x2)] = E 



E (^D 

+00 4 

3^ 



+ 2 



-rt 



-2 _^ 2x\ 



-3rt/2 



2x x"- 



+ 



zo(dx). 



Proof. Since r is constant, we have Ql = rl, where 1 is the constant mapping, which is equal 
to 1. Lemma [33] gives 

1 



nNt 



rE 



E fi^t) 

ueVt 



nfiYt)], 



where Y is generated by 



Afix) 



f{x) + 2r(f(- 



Finally, we complete the proof using [28, Theorem 4]. 



□ 



Now, let us talk about the rate of convergence. To estimate the distance between two random 
measures, we will use the Wasserstein distance B0]|44l : 

Definition 5.8 (Wasserstein distance). Let /xi and ^2 two finite measures on a Polish space 
{F, dp), the Wasserstein distance between and ^2 is defined by 



Wdpini,fJ.2) = inf 



FxF 



dF{xi,X2)U{dxi,dx2), 
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where the infimum runs over all the measures U on F x F with marginals /xi and fi2- In 
particular, if Hi and H2 are two probability measures, we have 

where the infimum runs over all two random variables Xi , X2, which are distributed according 
to pi,p,2- 

So, if Ml, M2 are two random measures then 

Wd{C{Mi),C{M2))=mi E[d{Mi,M2)], 

where the infimum is taken over all the couples of random variables (Mi, M2) such that Mi ~ 
£(Mi) and M2 ~ C{M2), and d is a distance on the measures space. Here, we consider 
d = VF|.|. It is the Wasserstein distance on {E, | • |). We have 

Theorem 5.9 (Quantitative bounds). If r is constant, then we have, for all t > 0, 




where Z^, are distributed as Z and start from 5.x and 5y. 

This result does not give a bound for Wiy,. I [C {Zt/E[Nt]) , C (W n)) or Ww^.^ (£ (Zt/iVt) , £ (vr)), 
where vr is the limit measure of Corollary 15.31 

proof of Theorem 1579] By homogeneity, we can see our branching measure Z as a process in- 
dexed by a Galton-Watson tree 0|. For our coupling, we take two processes indexed by the 
same tree. More precisely, as the branching time does not depend on the position, we can set the 
same times to our two processes. Let T = UngN{l; 2}" representing cells that have lived at a 
certain moment. Let {du)ueu be a family of i.i.d. exponential variables with mean 1/r, which 
model the lifetimes. We build Z^ and Z^ by induction. First, for all t £ [0, ^0), Xf = x + t and 
= y + t. We set a{i/}) = 0. Then, for all n G T and k G {1, 2}, we set a{uk) = a{u) + d„ 

and 

yt G [a{uk), a{uk) + Xf = + t - a{uk) 

and Yi"^ = Y^^^^^_/2 + t - a{uk). Finally we have Vt = {u e T \ a{u) < t < a{u) + du] 
and 

Zf = ^ 6xf and Z^ = ^ .5y.. 
u&Vt «eVt 
We observe that, for any cell u, the trajectories of X" and ai^e parallel (because they ai^e 
linear). When a branching occurs, X^neVt \-^t ~ ^t"l constant. Hence, we easily deduce that 

^ \Xf-Yl^\ = \x-y\. 
ueVt 
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Finally we have, for all t > 0, 

H^i.i(zf,zr)< j2 \^r-Yt^\ 

< \x — y\. 

Dividing by KlNt] = e"*"*, we obtain the first bound. For the second bound, a similar computa- 
tion gives 



Ww,, U Uf K^ kf 



< E 



1 



\x - y\ 



The process {Nt)t>o is know to be the Yule's process. It is geometrically distributed with pa- 
rameter e"*^'*, so we have 



E 



1 



rte' 



-rt 



1 



-rt 



It ends the proof. 



□ 



Remark 5.10 (Generalisation of Theorem 15. 9K In the proof of Theorem 1579] we only need that, 
for all n e N*, 6* G [0, l],t > 0, and x,y £ E 



Y,\F^'\Xt,e)-F^'\Yt,e)\<\x-y\ 
i=i 

where X, Y are generated by G and start respectively from x, y. For instance, we can consider 
that X is a continuous levy process and the division is a sub-critical fragmentation; namely 

Vx eE,yke N*,Vi < k, Fj''\x,e) = e'^x, 

where {&j)j,k a family of random variable verifying 



^ G)' < 1 and Vj G {1, . . . , k], O'^ G [0, 1]. 

Even if we do not find an explicit bound, we are able to prove a Wasserstein convergence. 
Lemma 5.11 (Wasserstein convergence). Under the assumptions ofTheorem \5.9\ we have 



lim Wy\ I t;-, vr 



in probability . 



Proof. As X 1-^ 1 + X is a Lyapounov function for the auxiliary process, we have 

, lim (/) = 7r(/) in probability, 

for all function / such that |/(x)| < C(l + x). The convergence also holds in distribution. By 
the Skorohod's Theorem, in another probability space, we have, 

lini ^(/) =7r(/)a.s. 
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for all continuous bounded function and for f{x) = x. This convergence is equivalent to the 
Wasserstein convergence. Thus, by a classical argument of discreteness (Varadarajan Theorem 
type), we get, 

lim W\.\ (^,it\ =Oa.s.. 

Hence, in our probability space we get that limf_s._|_oo W\.\ {Zt/Nt, vr) = in distribution. And 
as the convergence is deterministic, we get the result. □ 



5.3. Asymmetric mitosis : Macroscopic approximation. Now, we do not assume that the 
division is symmetric. We assume that 0) = F^^ {9)x and F2'^\x, 6) = {1 — F~^{6))x. 

We recall that F{x) = 1 — F(l — x). In this case. Equation ^ becomes 

dtn{t, x) + d^n{t, x) + r{x) n{t, x) = 2E[^r{x/Q)n{t, x/Q)], 

where n{t, .) is the density of X^. In particular, we deduce that the following P.D.E. has a weak 
solution: 

r+oa 

dtn{t,x) + dxn{t,x) + r{x)n{t,x) = / b{x,y)n{t,y)dy 

Jx 

where b verify the following properties: 



b{x, y) > 0, b{x, y) = for y < x 

r+oo 

b{x, y)dx = 2r{y) 



+ 00 

xb{x, y)dx = yr{y) 
b{x,y) = b{y - x,y). 

This equation was studied in |[36l . Here, 

b{x,y) = -r{y)g{-), 
y y 



(13) 
(14) 

(15) 
(16) 

(17) 



where g is the weak density of F. We easily prove the equivalence between to verify ([TT] ) and 
([13]- [HI). Our aim in this section is to describe the limit of the fluctuation process. It is defined 
by: 



Vt G [0,r],Vn G N*, 7?J"^ = 



Tlieorem 5.12 (Central limit Theorem for asymmetric size-structured population). Let T > 0. 

(n) 

Assume that tjq converges in distribution and that 



E 



sup / (l + x^) X[,"^(dx) 

n>l JE 



< +00. 



(18) 
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Then the sequence {ri^^^)n>i converges in ©([O, T], C~'^'^) to the unique solution of the evolu- 
tion equation: for all f G C'^'^, 

vtif) = m{f) (19) 

+ fix) + r{x) ^ f{qx) + /((I - q)x)F{dq) - f{x)j Vs{dx)ds 

+ Mt{f), 

where M(/) is a martingale and a Gaussian process with bracket: 

{M{f))t = / / 2f'{x)f{x) + 2r{x) / {f{qx) - f{x)fF{dq) Xs{dx)ds. 
Jo Jo Jo 

And C^'^ is the set of functions, such that f, f , f" vanish to zero when x tends to infinity. 
C~'^'^ is its dual space. 



Lemma 14.11 gives 



where for any / smooth enough, 

Vl")(/) = fj^^ nx)+T{x) i^J^' f{qx)+f{{l-q)x)F{dq)-f{x)^ r,t\dx)ds, 
and M^") is a martingale with bracket: 

(M(")(/))i = f r^2r{x) f\f{qx) - f{x)fF{dq) xW(dx) ds. (20) 
JO Jo Jo 

As the set of signed measure is not metrizable, we can not adapt the proof of Theorem 11.21 
Following 12911431 . we consider 7?(") as an operator in a Sobolev space, and use the Hilbertian 
properties of this space to prove tightness. See for instance 021 for condition to prove tightness 
on Hilbert spaces. Let us explain the Sobolev space that we will use. Let p > and j G N. The 
set W^'^ is the closure of , which is the set of functions of class C°° from into M with 
compact support, embedded with the following nomi: 

The set W^'^ is an Hilbert space and we denote by W~^'^ its dual space. Let C^'^ be the space 
of function / of class such that: 

f^^Mx) 

Wk < j, lim ^ = 0. 

a::->+oo 1 -\- xP 

We embed it with the following norm: 

v/GC^'^ ll/lb. = E-PTv^- 

The set C^'^ is also a Banach space and we denote by C~^'^ its dual space. These spaces verify 
the following continuous injection |[29l [D : 

Ci,P C W^'P^^ and VF^+^'P c C^'^. (21) 
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Or equivalently, if / is smooth enough, 

< C\\f\\c,„ and WfWc.p < C\\f\\wi+i„. 

The first embedding/inequality prove that the tightness in W^''^'^^ implies the tighness in C^'^. 
The second is useful for some upper bounds. For instance, we have 

Lemma 5.13. If {ek)k>i is a basis ofW^'^ then we have, for all k > and x ^ E, 

Y.eu{xf<C{l + x^). 

k>l 

Proof. Sx : f ^ f{x) is an operator on W^'^. We have, for all / G W"^'^, 

%f\ < (1 + ^)||/||c<u < C{1 + < C{1 + 

But, by Parseval's identity we get. 



k>l 



which completes the proof. 

We introduce the trace (((M(")))j)^^^ of (m|"^)^^^. It is defined such that 



Mi")||2,_.,.-((MW)) 



in)\ 



t>0 



is a local martingale. Then since 



M 



(n) 



w- 



k>l 



□ 



where {ek)k>i is a basis of W'^'^. Then by (l20l ). we get 

, , ft /'+00 fl 

((M(")))t = ^/ / 2r{x) {eu{qx)-ek{x)fF{dq)l^^p\dx)ds. 
g^Jo Jo Jo 

Now, we first prove the tightness of (?7^"^)n>i then Theorem 15. 121 
Lemma 5.14. (r/")„>i is tight in B{[0, T], W~^''^) 

Proof By |[24l Theorem 2.2.2] and Il24l Theorem 2.3.2] (see also ||29] Lemma C]), it is enough 
to prove 

(1) e[sup,<J|77^||2^,„2,i] < +00, 

](5 > such that for each stopping times Sn bounded by T 



(2) Vn G N, Ve, p > 0, 



n-s>+oo 0<u<5 



lim sup sup I 

n-s-+oo 0<«<5 



Sn+U 
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For the first point, using lemma ISTTl there exists C > such that 



k>l k>l k>l 



<CxS"^(l + x) 



Then, since 



M 



in) 



w- 



k>l 



Doob's inequahty and (ITSl) gives 



E 



sup 

te[o,t] 



M 



(n) 



w- 



<C', 



where C > 0. Then there exits C" > such that 



m 

And as 



(n) 



VK-2.1 



< 



(n) 



VK-2.1 



V 



(n) 



the Gronwall Lemma gives 



+ 



w- 



E 



V 



(n) 



+ 



M 



(n) 



VK-2.1 



<C" + 



V 



(n) 



< C / sup 



(n) 



w- 



ds, 

!,1 



sup 



(n) 







rS„+u 




2 




< E 


K' f sup 

JSu s<T 





for a certain constant K. Finally for the second point, we have 

'Jn~\~U On 

< K"u. 

Here K' , K" are two constants. Using the Markov-Chebyshev inequality, we prove the Aldous 
condition. We similarly prove that ((M^"))) verifies the Aldous condition. We deduce that 
(?7("))„>i is tight. □ 



Proof of Theorem \5.12\ Let M be a continuous Gaussian process with quadratic variation veri- 
fying, for every / G C^'O (c VF^'^) and t G [0, T], 

{Mif))t = J2 / 2r(x) / ifiqx)-fix)fFidq)Xsidx). 
p^iJo Jo Jo 

Since there exists C f such that 



C 



V/gC^'O, sup |m(")(/)|<-4, 

te[o,T] V"- 

and (M("))t converges m law to {M)t, then by 1231 Theorem 3.11 p.473], M(")(/) converges 
to M(/) in distribution, as n tends to oo. 
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By Lemma [5.14l and (|2TI ) , the sequence (?7^"^)n>i is also tight in C Let r] be an accumula- 
tion point. Since its martingale part M in its Doob's decomposition is almost surely continuous, 
then rj is also almost surely continuous. Hence, is a solution of ([T9l ). Using Gronwall's in- 
equality, we obtain the uniqueness of this equation, in C([0, T], C~^'*^), up to a Gaussian white 
noise M. We deduce the announced result. □ 

6. Another two examples 

6.1. Space-structured population model. Here, we study an example which can models the 
cells localisation. One cell moves following a diffusion on C W^, d>l, and when it dies, its 
offspring is localised at the same place. Hence, in all this section the branching is local; that is 

VA; > 0,Vi < fc,Vx G E^'iO G [0,1], F^'"\x,6) = x. 

6.1.1. Branching Omstein Uhlenbeck. In this subsection, we consider the model of |[T5l Exam- 
ple 10]. Assume that G is given by 

G/(x) = Va/(x)-5x.V/(x), 

where o", (7 > 0, f is smooth, x G W^. Also assume that the division is dyadic, that is p2 = 1, 
with rate 

r{x) = hx^ + a, 

where a, 6 > and a or 6 is not null. Here x"^ = = x.x. If 5 > \/26 then we add the 
following notations: 



9-^9'- 2ba' 



and a = \ g"^ — 2ba'^ . 



2a2 

We also denote by tToo the Gaussian measure whose density is defined by 



a \ fa 



From our main theorem, we deduce 

Corollary 6.1 (Limit theorem for an branching Omstein-Uhlenbeck process). If g > a\/2b and 
Xq = X G M'^ then, for any continuous and bounded f, we have 

t^+oc Nt /igd y TToo {dy) 

in probability. In particular, 

^here A = ^"^g^-^"'-^ + 

a is the Malthus parameter 

Proof. \fV : X ^ e^^^ then it is an eigenvector of Q, which is defined for any smooth / by 

gf{x) = Gf{x)+r{x)f{x). 
We conclude the proof using Theorem 1 1.1 1 and Lemma 1331 □ 
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Remark 6.2 (Another eigenelement). Note that ifV2 x ^ e^^^'^ then it is an eigenvector ofQ, 
associated to the eigenvalue 

A2 = + a. 

But in this case, the auxiliary process is not ergodic and we are not able to deduce any conver- 
gence from our main theorem. 

6.1.2. General case. Let us assume that G is the generator of a diffusive Markov process. If 
the state space E is bounded then we can find sufficient conditions to the eigenproblem in 
section 3] and f38l Theorem 5.5]. For instance, under some assumptions, we have 



Ao = lim In supE 



Nt I Xo" 



If E is not bounded then we can see 112111411 . This example is developed in fT5l . They prove a 
strong law of large number, which is close to Theorem 11.11 

6.2. Self-similar fragmentation. Self-similar mass fragmentation processes are characterised 
by 

• the index of self-similarity a G M; 

• a so-called dislocation measure v on S = {s = (sj)jgN I limj_j.+oo Si = 0, 1 > Sj > 
Si > 0, Vj < i} which satisfies 



1^(1,0,0, ..) = and j {I - s)u{ds) < +oo. 



If v{S) < -\-oo then the dynamics is as follows: 

• a block of mass x remains unchanged for exponential periods of time with parameter 

x"z^(5); 

• a block of mass x dislocates into a mass partition xs, where s G 5, at rate u(ds); 

• there are finitely many dislocations over any finite time horizon. 

The last point is not verified when iy{S) = +oo. In this case, there is a countably infinite 
number of dislocations over any finite time horizon. So, when < +oo, our setting capture 
this model with the following parameters: 

G = 0, r(x) = xM5), f'Y.Pkix)T.fi^J^^^)de= f Y.fis.xf^'^''^ 

•^0 k>o j=i ■^■^ i>0 



Hence, in this case we have 



and V : X xP is an eigenevector. See ||8l for further details. 



i>0 \ i>0 
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